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We consider supersymmetric models with right-handed neutrinos where neutrino masses are purely 
Dirac-type. In this model, right-handed sneutrino can be the lightest supersymmetric particle and 
can be a viable candidate of cold dark matter of the universe. Right-handed sneutrinos are never 
thermalized in the early universe because of weakness of Yukawa interaction, but are effectively 
produced by decays of various superparticles. We show that the present mass density of right- 
handed sneutrino can be consistent with the observed dark matter density. 

PACS numbers: 14.60.Pq, 12.60. Jv 98.80.Cq, 95.35.-|-d 


In recent years, various experiments have confirmed 
the phenomenon of neutrino oscillation. (See, for ex¬ 
ample, Hiiii .) Those results strongly suggest 
very small but non-vanishing neutrino masses. This fact 
raises serious problems because the non-vanishing neu¬ 
trino mass is not allowed in the standard model of parti¬ 
cle physics and also because suggested values of neutrino 
masses are extremely small. The easiest way of generat¬ 
ing neutrino masses is to introduce right-handed neutri¬ 
nos; with this extension, Yukawa couplings of neutrinos 
may exist. Consequently, neutrinos can acquire masses 
after electroweak symmetry breaking. 

Even with right-handed neutrinos, there are two dif¬ 
ferent classes of scenarios for generating neutrino masses. 
Probably, more popular one is with Majorana masses for 
^ht-handed neutrinos, i.e., so-called “seesaw” scenario 
|6|. In this scenario, smallness of the neutrino masses is 
explained by the Majorana masses of right-handed neu¬ 
trinos which are much larger than the electroweak scale. 

Small neutrino masses can be, however, realized with¬ 
out seesaw mechanism. With vanishing Majorana masses 
of the right handed neutrinos, which may be the con¬ 
sequence of exact lepton-number symmetry, neutrino 
masses become Dirac-type. As we will see, in this case, 
Yukawa coupling constants for neutrinos are, roughly 
speaking, 0(10“^^) or smaller to make the neutrino 
masses to be consistent with the results of neutrino- 
oscillation experiments. One might think that such small 
Yukawa coupling constants are unnatural. It is, how¬ 
ever, natural in ’t Hooft’s sense 0 since some symme¬ 
try (i.e., chiral symmetry in neutrino sector) is restored 
in the limit of vanishing neutrino Yukawa coupling con¬ 
stants. In this model, however, right-handed neutrinos 
are mostly irrelevant for collider experiments and cos¬ 
mology since their interaction is extremely weak. 

In supersymmetric models, which is strongly motivated 
as a solution to various problems in the standard model of 
particle physics, like hierarchy and naturalness problems, 
situation changes. In particular, superpartners of right- 
handed neutrinos may play an important role in cosmol¬ 
ogy. When small neutrino masses are purely Dirac-type, 
masses of right-handed sneutrinos are dominantly from 
effects of supersymmetry (SUSY) breaking. Then, one 
should note that the lightest superparticle (ESP) may 


be the (lightest) right-handed neutrino vr. Importantly, 
since the LSP vr becomes stable by the i?-parity con¬ 
servation and also is very weakly interacting, it can be 
a viable candidate of cold dark matter (CDM) provided 
that its relic density is the right amount. 

In this letter, we consider the minimal supersymmet¬ 
ric standard model (MSSM) with three generations of 
right-handed (s)neutrinos where small neutrino masses 
are purely Dirac-type. In particular, we study the case 
where the LSP is the lightest right-handed sneutrino and 
see if the relic density of vr can become consistent with 
the present CDM density. Since interaction of right- 
handed sneutrino is very weak, it is not thermalized in 
thermal bath. Even in this case, some of decay (and scat¬ 
tering) processes produce vr. As we will see, in some 
parameter region, the density parameter of right-handed 
sneutrino can be 0(0.1), which is consistent with the 
CDM density suggested by the WMAP 0 : 

Q 7,2 — 0 ('ll 


where /iioo is the Hubble constant in units of 
100 km/sec/Mpc. 

Let us first introduce interaction and mass terms in 
Lagrangian. The important part of superpotential is 

W = VrHuLu’r + mHuHd, (2) 

where Hu = {Hu,H^) and Hd = {Hd,H^) are up- and 
down-type Higgses, and L = {Dr, If) the left-handed lep¬ 
ton. (In this letter, “hat” is for superfields while “tilde” 
is for superparticles with odd i?-parity.) Here, we omit 
flavor indices for simplicity. With this superpotential, 
neutrino mass is generated after electroweak symmetry 
breaking: rriu = yv{Hu) = yvVS\a.f3, where v ~ 174 GeV 
and tan/3 = (i7°)/(iJ°). Thus, the Yukawa coupling con¬ 
stant yu is determined once the neutrino mass is fixed: 


yu sin/3 = 3.0 x 10 x 


2.8 X 10-3 eV2 


1/2 


( 3 ) 


For simplicity, we consider the case where tan /3 is rela¬ 
tively large. In this case Hu behaves like the standard- 
model Higgs; Hu — Hsm- The lightest Higgs boson 
h is contained in Hsm, and we take its mass to be 



2 


rrih = 115 GeV in our study. On the contrary, Hd plays 
the role of heavy Higgs doublet. In addition, we ap¬ 
proximate mass eigenstates of the superparticles by left- 
and right-handed sleptons, gauginos, and Higgsinos. Ef¬ 
fects of left-right mixing are taken into account by mass- 
insertion method. The mass terms for superpaticles are 
given by 


1 

2 1 


rrigBB + niy^WW + ^jLnHuHd + h.c. 
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where B and W are gauginos for U{1)y and SU{2)r 
gauge groups, respectively. Notice that all the terms ex¬ 
cept for Higgsino mass term are from SUSY breaking 
and that m- and to? are identical up to Z?-term contri- 
bution. Furthermore, we introduce soft-SUSY breaking 
tri-linear coupling which will become important in the 
following discussion 

C-A = h.c. (5) 

For later convenience, we parameterize 

+ = a^y^rni, (6) 


Decay rates for processes which become important for 
the calculation of are given by 


^HO^Pri?l ~ ^H+^pRlf — 
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where we have neglected lepton masses. Here, gi and 52 
are gauge coupling constants for the U{1)y and SU{2)l 
gauge groups, respectively, and, for the process x vrij^ 
/3f is given by 


where represents the D-term contribution. In ad- 

dition, we adopt the GUT relation among the gaugino 
masses: to^ ~ O.Stoj^. 

Now let us discuss the production of the LSP right- 
handed sneutrino vr via various decay processes. Since 
we found that the production via scatterings is subdom¬ 
inant and can be neglected, we consider here only the 
decay processes x —> VRy. We assume that particles 
X and y are in chemical equilibrium. Their distribu¬ 
tion functions at the temperature T, then, are given by 
f{E) = (e^/^ ± 1)“^, where the positive and negative 
signs are for fermions and bosons, respectively. 

The evolution of the number density of vr is governed 
by the Boltzmann equation 


P( = ‘^{rnlR + ml)ml + (to?^ - to^)^], (16) 

with rux and ruy being the masses of the particles x and 
y, respectively. In addition, vt is temperature-dependent 
Higgs VEV. We approximate the Higgs potential in ther¬ 
mal bath as 0 

Ut ^ 

+ g;j^(2TO^-|-TO^ + 2TO()T^|i?SMp, (17) 

where T is cosmic temperature. Then we minimize Vt 
to obtain vt — {Hu)t- 

In solving Eq. 0 it is useful to define the yield variable 


BpR — Gdecay; (7) 

where “dot” represents time-derivative and H is the Hub¬ 
ble parameter. The source term from the decay process 
is written in the form 

■ 5 -^ f dPk 

Gdecay = J + i)Tx^pRyfx{i E fy)k^ ,{&) 

where 7 ^ = ruxt Lorentz factor and 

(25^: + 1) is the spin multiplicity of x. Furthermore, 
(1 ± fy)k^ is the averaged final-state multiplicity factor 
for fixed value of initial-state momentum. (Here, the 
positive and negative signs are for bosons and fermions, 
respectively.) 


( 18 ) 

where s = is the total entropy density of the 

universe. (In this analysis, we take the effective number 
of massless degrees of freedom as g^ = 106.75, since the 
production of vr becomes effective at the temperature 
lower than most of superparticle masses.) 

Using the relation T = —HT, Yp^ is given by 

p'^ma.x 

Yp^{T) = J^ -^dT. (19) 

Once the yield variable is obtained, we can also calculate 
the mass density of sneutrino using the relation PprIs = 
mp^Yp^. Gomparing with the present value of the critical 
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density pcrit, which is given by [/9crit/s]now 3.6 x 
10“® GeV, we obtain density parameter of right-handed 
sneutrino = {pvr + Pt* )/pcrit- 

Importantly the present value of Yor is insensitive to 
the maximal temperature T„iax of the universe. To see 
this, it is instructive to roughly estimate Yor . Neglecting 
the Lorentz factor, let us approximate the decay term as 
Gdecay "^^^J/^msuSY^^suSY, where msuSY and nsuSY 
are typical mass scale and number density of (parent) 
superparticles, respectively, and A^mode is the number 
of possible decay mode. The number density nsuSY is 
^ when the temperature is higher than msuSY) and 
is exponentially suppressed when T ^ msuSY- The in¬ 
tegration in Eq. m is then dominated at the tempera¬ 
ture of T ~ TOsuSY- Consequently, we obtain Yp^(T ^ 

wsusy) ~ ^ ^ 2.4 X10^® GeV being 

167rg,' msuSY 

the reduced Planck scale. As one can see, the present Yp^ 
is insensitive to thermal history for T ^ tosusy- Thus, 
does not depend on, for example, reheating temper¬ 
ature after inflation. We can also estimate the density 
parameter; using the above estimation of we obtain 
^ORhl^^ - O(ioy3) X With 

this naive estimation, ^Or becomes smaller than the cur¬ 
rently observed CDM density. In some case, however, 
^Pr becomes much larger, as we see below. 

Now, we are at the position to quantitatively estimate 
^Pr. We have numerically evaluated the yield variable 
using Eq. m for several choices of parameters. The relic 
density of right-handed sneutrino strongly depends on 
the neutrino Yukawa coupling constant which is related 
to the neutrino mass. Neutrino-oscillation experiments 
determine only the mass-squared differences of neutrinos; 
here we adopt ili 

Maton. = 2.8xlO-®eV^ (20) 

Msolar = 7.9xlO-®eV^ (21) 

Eirst, let us consider the case where the masses of the 

neutrinos are hierarchical. In this case, [AmJ] al¬ 

most corresponds to the mass-squared of the heaviest 
neutrino which we call third generation neutrino. In this 
case, the largest Yukawa coupling constant is given by 
~ 3 X 10“^®, and other Yukawa coupling constants 
are much smaller. (Here and hereafter, the superscript 
“(*)” indicates that is for neutrino in z-th generation.) 
As a result, production of £'_r is dominated by processes 
where the third-generation (s)neutrino is related. 

As show in in Eqs. 0 - itT^ . there exist various 
decay processes which produce vr. Among them, Hig- 
gsino decay process dominates the vr production when 
the effects of the tri-linear scalar couplings are negligi¬ 
ble. In this case, however, Q.p^ becomes too small to 
be consistent with observation of the GDM relic density. 
Indeed, when = 100 GeV, for example, we found 
VLp^ = 0.004 - 0.001 for /r// = 200 GeV - 1 TeV, which 
is much smaller than the present CDM density. 

If some of other decay processes are effective, VLp^ can 
become significantly larger. In particular, when the tri- 
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FIG. 1: as a function of for a^, = 1 — 6. Here, we 

take mpR = 100 GeV, = 300 GeV, and rh = 150 GeV. 
The shaded band corresponds to the CDM density suggested 
by the WMAP. 


linear coupling constant Ai, is non-vanishing, decays of 
various superparticles produces vr sufficiently. 

One possible enhancement can be due to the process 
W —> Vr + ■ ■ ■■, if the mass of left-handed sneutrino 
Vr is close to mpj^, production of vr is enhanced be¬ 
cause the left-right mixing of sneutrinos becomes larger 
as becomes smaller. (See Eqs. iHTHl — llTKll .) 

In fact, Q,p^ can be of 0(0.1) with a mild degeneracy of 
VR and VL- To see this, in Eig.^] we show as a 

function of mp^. We can see that 10 — 20 % degeneracy 
is enough to realize fi/i-CDM even if a,y < 3. Here, we 
take = 300 GeV. If we increase the Wino mass, ^p^ 
decreases since the temperature-dependent VEV vr be¬ 
comes suppressed at high temperature. We also note that 
the process —> vrt'^ gives additional contribution to 

^Pr, which may become sizable when tan/3 is large. We 
found that ^Ip^ can be ~ 40% larger for tan/3 = 55. 

Even without mass degeneracy between vr and vr, vr- 
GDM is realized if ^ 1. Although such a very large 
value of an may not be realized in simple supergravity 
models, it is phenomenologically viable. 

Next, we turn to consider the case of degenerate neu¬ 
trino masses, where there is another possibility of en¬ 
hancing O.Pr. In this case, the relation y^^ ~ y^'^ ~ y^^ 
holds and all three generations of right-handed sneutri¬ 
nos may be effectively produced. The point is that the 
neutrino Yukawa coupling constants can be much larger 
than ~ 3 X 10“^®, and hence Q.p^ can be more enhanced 
than the hierarchical case. />/j-CDM can be realized if the 
neutrino masses are 0(0.1 eV) even if there is no other 
enhancement. 

So far, we have only discussed the production of vr 
by the decay of superparticles in chemical equilibrium. 
In general, vr can be also produced after superparti- 
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cles freeze out. The lightest superparticle in the MSSM 
sector, which we call MSSM-LSP, is very long-lived be¬ 
cause of the weakness of the neutrino Yukawa interac¬ 
tion. Thus, decay of the MSSM-LSP occurs after its 
number is frozen, which gives another source of vr. 
The total density parameter is then given by = 

+ mMssM ^MSSM, where tumssm is the MSSM-LSP 
mass. Here, the first term is evaluated with Eq. m 
while the second one is from the decay after the freeze 
out; Hmssm is the expected density parameter of the 
MSSM-LSP for the case where the MSSM-LSP is sta¬ 
ble. Notice that Hmssm is model-dependent and may be 
smaller than the WMAP value; if Hmssm is large, the 
parameter space to realize izfl;-CDM is enlarged. 

Here we discuss the effects of the MSSM-LSP decay 
on big-bang nucleosynthesis (BBN). If the MSSM-LSP 
decays after BBN starts, it may cause hadro- and photo¬ 
dissociation of light elements, which may spoil the success 
of standard BBN scenario. In order to realize i/^j-CDM 
with enhanced left-right mixing, is required to be rel¬ 
atively light and then it may be the MSSM-LSP. If so, we 
found that, for = 100 GeV, its lifetime becomes so 
long that light-element abundances become inconsistent 
with observations due to hadro-dissociation processes m 
unless rrioL ^ 140 GeV (with Oj, > 6). On the other hand, 
very degenerate sneutrinos (rnoj^ — ^ a few GeV) is 

another possibility to avoid the difficulty by suppressing 
the energy release from i/r decay. Otherwise, the BBN 
constraints can be also avoided if the MSSM-LSP is the 
lightest neutralino, or if i?-parity in the MSSM sector is 
(very weakly) broken. 

Finally, we comment on the decay of heavier right- 
handed sneutrinos (denoted by ^r)- When the mass 
difference between and vr is large enough, may 
decay, for example, into H^L', and the produced Hig- 
gsino decays subsequently as — > vrL. In this case, 

lifetime of the long-lived particles can be ~ 1 sec or 
shorter; such decay processes are cosmologically safe m- 
On the contrary, if decays only through the pro¬ 
cesses of —> i>R + ■ ■ ■ due to kinematical reason, 

the lifetime becomes much longer than the present age 
of the universe. For instance, when mp^ = 100 GeV, 
mp'^ = 150 GeV, and = 200, 300, and 400 GeV, the 


lifetime of is estimated as Tp'^ ~ 5 x 10^^, 5 x 10^®, 
and 2 x 10^® yr, respectively. (Here, we approximated 
that all the neutrino Yukawa coupling constants are 
10“^^, and neglected tri-linear scalar couplings.) Even 
in this case, some (tiny) amount of heavier right-handed 
sneutrinos decay into charged particles until today. As¬ 
suming the density parameter of i>' to be 0(0.01), flux 
of the (primary) charged particles is estimated to be 
T ~ 10-16 cm-^sec-istr-iGeV-i x ■ The 

emitted charged particles scatter off the background ra¬ 
diation and produce energetic photons. However, we 
expect that such small flux is consistent with the ob¬ 
servations since the above flux is much smaller than 
the observed y-ray flux with the energy 10 — 100 GeV: 
yr(obs) ^ xo“8 — 10-16 cm-^sec-istr-iGeV-i [I^ . 

In summary, in this letter, we have considered super- 
symmetric models with the right-handed-sneutrino LSP, 
where masses of neutrinos are purely Dirac-type. We 
have seen that, in such models, ^p^ can be of 0(0.1) and 
right-handed sneutrino can be CDM. In particular, if vr- 
CDM is realized in the parameter region where vr and vr 
are quite degenerate, left-handed sleptons become fairly 
light. Such light sletons are interesting targets of the fu¬ 
ture collider experiments. On the contrary, if i/ij-GDM 
is realized with degenerate neutrino masses, the (active) 
neutrino mass should become ~ 0(0.1 eV). Such a sce¬ 
nario may be tested by using matter power spectrum 
from Ly-a forest data jl^ • 

If Vr is the LSP, there are much more to be discussed 
from phenomenological and cosmological points of view. 
Even though the LSP is vr, the MSSM-LSP looks like 
stable particle in collider experiments. However, in this 
case, the MSSM-LSP may be charged or even colored. 
In addition, cosmology may be also significantly affected 
if the right-handed sneutrino is the LSP; one of the ex¬ 
amples has been discussed in this letter. We also note 
that, if the initial amplitude of vr is ~ 10® GeV, coher¬ 
ent mode of the right-handed sneutrino may become the 
CDM. We believe that the scenario with right-handed- 
sneutrino LSP provides new and rich phenomenology. 
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